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Abstract 

We present a stochastic analysis of a data set consisiting of 10 6 
quotes of the US Doller - German Mark exchange rate. Evidence is 
given that the price changes x(t) upon different delay times r can be 
described as a Markov process evolving in r. Thus, the r-dependence 
of the probability density function (pdf) p(x, r) on the delay time r 
can be described by a Fokker-Planck equation, a gerneralized diffusion 
equation for p(x,r). This equation is completely determined by two 
coefficients Di(x,t) and D2(x,t) (drift- and diffusion coefficient, re- 
spectively). We demonstrate how these coefficients can be estimated 
directly from the data without using any assumptions or models for 
the underlying stochastic process. Furthermore, it is shown that the 
solutions of the resulting Fokker-Planck equation describe the empir- 
ical pdfs correctly, including the pronounced tails. 
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1 Introduction 

Since L. Bachelier's pioneering work dating back to 1900 [[|], the complex 
statistical properties of economic systems have attracted the attention of 
many researchers and an extensive literature has evolved for modeling fluc- 
tuations of financial markets. Traditionally, fluctuations of financial assets 
were viewed and modeled as random variables. Well known examples are the 
ARCH-type models (see for example, p], ^ [§]) and stochastic volatility (SV) 
models (|, 

Since advances in computer technology have made high frequency data 
available, many physicists have joined the field adapting methods from sta- 
tistical physics. One line of studies within econophysics focusses on the sta- 
tistical properties of financial time series such as stock prices, stock market 
indices or currency exchange rates. Rather than comparing the predictions of 
models with the various aspects of emprical data (which is the traditional ap- 
proach), physicists try to extract information about the stochastic processes 
governing financial markets from an analysis of emprical data. An overview 



of recent developments in this field can be found in references |9j , |TI| and 



TTj] . In spite of all the efforts spent, some of the most basic questions con- 
cerning the statistics of financial assets remain unsolved. In particular, the 
mechanism leading to the fat-tailed (leptokurtic) probability distributions of 
fluctuations on small time scales is still unknown. Compared to a Gaussian, 
the pdfs of those fluctuations express an unexpected high probability for large 
fluctuations (see figure |13|). Quantifying the risks of such large fluctuations 
is of utmost importance for the risk management and the pricing of options. 

Fluctuations of financial time series y(t) are usually measured by means of 
returns, log-returns or, equivalently, by price increments. Here, we consider 
the statistics of the price increment x(t) over a certain time scale r, which 
is defined as: 

x(t) = y(t + T)-y(t). (1) 

We suppressed dependence of the price increment x(r) on the time t since 
we assume the underlying stochastic process to be stationary. 

In a recent letter ||12|| , we demonstrated that the mathematics of stochas- 
tic processes is a useful tool for empirical investigations of the time scale 
dependence of the proabability density function (pdf ) p(x, r) of the price 
increment x on the time scale tQ. It was shown how the equations governing 

1 The correct notation for the probability density function of the stochastic variable 
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the underlying stochastic process can be extracted directly from the empir- 
ical data, provided that the price increment x obeys a Markov process. In 
particular, it is possible to derive a partial differential equation, the Fokker- 
Planck equation, which describes the evolution of the pdf p(x, r) in the scale 
variable r. Hence, the mathematics of Markov processes yields a complete 
description of the stochastic process underlying the evolution of the pdfs 
from Gaussian distributions at large scales r to the leptokurtic pdfs at small 
scales. 

Here, we extend our previous analysis presented in considerably. We 
show how the existence of a Markov process can be checked empirically and 
how the Fokker-Planck equation can be calculated directly from the data. 
Furthermore, we show how multiplicative and additive noise sources interact 
leading to the leptokurtic distributions and volatility clustering of the ex- 
change rates. Finally, our results are interpreted in the context of currently 
discussed features of financial time series. Since our method is not based on 
any models, we gain new insights into the mechanisms governing the statis- 
tics of economic systmes. In particular, we find new evidence for cascading 
processes in financial markets. 

The paper is organized as follows. Section ^| is devoted to a brief summary 
of the most important notions and theorems on Markov processes and their 
application to the analysis of empirical data, while section ^| contains the 
main results of our analysis. We shall consider in detail a data set containing 
10 6 samples of the DM / US-Dollar exchange rate from the one-year period 
October '92 to September '93. An interpretation of our results is finally given 
in section 0. 



2 Mathematics of Markov Processes 

The following section briefly summarizes the notions and theorems which will 
be of importance for our statistical analysis. For further details on Markov 
processes we refer the reader to the references JR] and [IS|. 

Fundamental quantities related to Markov processes are conditional prob- 
ability density functions. Given the joint probability density p(xi, T\\ x 2 , t 2 ) 
for finding the price increments x\ ( := x{t\) ) at a scale t\ and x<i at a scale 

x(t) is p(x(r)). We write p(x, r) in order to indicate that the pdf of x(t) depends on the 
parameter r. 
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r 2 with 7~i < r 2 , the conditional pdf p(x\, Ti\x 2 , r 2 ) is defined as: 

p{Xl,T 1 \X2,T 2 ) = r . (2) 

P{X2,T 2 ) 

p(x\, Ti\x2, t 2 ) denotes the conditional probability density for the price incre- 
ment x\ at a scale T\ given an increment x 2 at a scale r 2 . It should be noted 
that the small scale T\ lies within the larger scale r 2 (see equation (|l|): the 
increments x\ and x 2 are calculated for the same time t). 

Higher order conditional probability densities can be defined in an anal- 
ogous way: 

/ i \ p(x 1 ,t 1 ;x 2 ,t 2 ;...;x n ,t n ) 

P(x 1 ,t 1 \x 2 ,t 2 ;...;x n ,t n ) = r . (3) 

P[x 2 ,t 2 ; ...;x N ,r N ) 

Again, the smaller scales T{ are nested inside the larger scales r i+1 (with the 
common reference point t). 

The stochastic process in r is a Markov process, if the conditional prob- 
ability densities fulfill the following relations: 

p(x 1 ,T 1 \x 2 ,T 2l X 3 ,T 3 ] . . .X N ,T N ) = p(xx,Tl\x 2 ,T 2 ) (4) 

with t\ < r 2 < r 3 < . . . < Tjy. As a consequence of (|j), each N-point prob- 
ability density p(x±, T\\ x 2 , t 2 ; . . . ; xn, Tat) can be determined as a product of 
conditional probability density functions: 

p(x 1 ,T 1 ;...X N ,T N ) = p(x 1 ,T 1 \x 2 ,T 2 )....p(x N _ 1 ,T N _ 1 \x N ,T N )p( y X N ,T N ) . (5) 

Equation (]5|) indicates the importance of the conditional pdf for Markov 
processes. Knowledge of p(x, t\xq, r ) (for arbitrary scales r and r with 
t < To) is sufficient to generate the entire statistics of the price increment 
encoded in the N-point probability density p(x\, T\\ x 2 , r 2 ; x 3 , r 3 ; . . . ; xn, Tn). 

For Markov processes the conditional probability density fulfills a master 
equation which can be put into the form of a Kramers- Moyal expansion 0: 

d 00 / d \ k 

- r —p(x,t\x ,t ) = 22 y-T^; J D k (x,T)p(x,T\x ,T ). (6) 

2 Note that, in contrast to the usual definition as, for example, given in [14] , we mul- 
tiplied both sides of the Kramers-Moyal expansion by r (see also equation (pf)). This is 
equivalent to the logarithmic length scale —ln( At / 40960s) used in fll^ . The negative sign 
of the left side of equation (||) is due to the direction of the cascade torward smaller scales 

T. 
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The Kramers-Moyal coefficients Dk(x, r) are defined as the limit At — > of 
the conditional moments Mk(x, r, Ar): 

£> fc (x,r) = lim M fe (x,r, At) , (7) 

At— >0 

+oo 

M k (x,T,Ar) = rf^~ / (5 - P (5, ^ - Ar|x, r) rfx. (8) 

— oo 

For a general stochastic process, all Kramers-Moyal coefficients are different 
from zero. According to Pawula's theorem, however, the Kramers-Moyal ex- 
pansion stops after the second term, provided that the fourth order coefficient 
D±(x,t) vanishes. In that case, the Kramers-Moyal expansion reduces to a 
Fokker-Planck equation (also known as the backwards or second Kolmogorov 
equation) : 

d \ d d 2 ) 

- t—p(x,t\x ,t ) = I -—D x {x,t) + —D 2 (x,t) *>p(x,t\x ,t ). (9) 

Di is denoted as drift term, D 2 as diffusion term. The probability density 
p(x, t) has to obey the same equation: 

- t^ P {x,t) = | -^D 1 (x,t) + ^D 2 (x,r) j p{x,r). (10) 

The Fokker-Planck equation describes the probability density function of a 
stochastic process generated by the Langevin equation (we use Ito's defini- 
tion) 

" r JW) = ^ r ) + y/D 2 {x,r)f(r) (11) 

where /(r) is a Langevin force, i.e. 5-correlated white noise with a gaussian 
distribution. Here, the increment x(t) at a fixed time t is generated by 
a stochastic process with respect to the continuous variable r. This kind 
of process expresses nothing else but a hierarchical, cascade-like structure 



connecting price increments on different time scales |L6 



In the case that the random force f(r) ist not Gaussian distributed, the 
Kramers-Moyal expansion does not reduce to the Fokker-Planck equation. 
From the more general Kramers-Moyal expansion (|6]), which is also valid 
for the probability density p(x,r), differential equations for the n-th order 
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moments can be derived. By multiplication with x n and integration with 
respect to x we obtain: 

a oo + °° / ft \ k 

-r— <x b (t)> = 2(-l)*yx B ^J D*(x,T)p(x,r)dt; 

Zj ( n _ fe ); y %n Dk{x,T)p{x,r)dv . 



k= 

n 



EAn(^% r )) • ( 12 ) 



3 Data Analysis 

The hypothesis concerning the Markovian properties of exchange rate data 
immediately fixes a framework for the analysis of the data. First, one has 
to give evidence of the Markovian properties according to equation Sec- 
ondly, the evolution of conditional probability densities p(xi, T\\x2, t 2 ) has to 
be specified on the basis of the stochastic evolution equation, eq. (|(J). To this 
end, we have to determine the conditional moments M^(x, r, At) at different 
scales r for various values of At. 

Practically, it is only possible to evaluate the lowest order coefficients. 
Therefore, we shall restrict our analysis to the coefficients of order one, two, 
and four. Approximating the limit At — > 0, we obtain the Kramers-Moyal co- 
efficients Dk(x, t). If the fourth order coefficient vanishes, the evolution equa- 
tion (|5p takes the form of a Fokker-Planck equation, according to Pawula's 
theorem. Since we consider a finite data set, we are not able to prove rigor- 
ously (in a mathematical sense) that D4 is zero but can only give hints for 
the validity of the assumption that the conditional probability density obeys 
a Fokker-Planck equation. 

In order to verify this assumption and our results for the Kramers-Moyal 
coefficients, we compare the numerical solutions of the Fokker-Planck equa- 
tions (^) and ( PD with the probability density functions determined directly 
from the data. 

For the rest of the article, prices and price increments are given in units 
of the standard deviation a of y(t). For the data set under consideration, a 
is 0.064DM. 
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3.1 Data Processing 

From the frequency spectrum (figure |l|) of the exchange rate y(t), it becomes 
evident that the data are dominated by white noise for high frequencies 
(above ~ 10 -2 Hz), comparable to the noise affecting any physical measure- 
ment. The original signal y(t) seems to be the sum of the underlying "real" 
signal s(t) and some additive white noise n(t): 



y (t) = s{t)+n(t). 



(13) 



The additive noise acts on the increments in a similar way (see equation 
([]])). Note that this noise differs from the dynamical noise of the Langevin 
equation (0). n(t) can be regarded as some randomness which is added 
seperately to the dynamical stochastic cascade process. In physics this kind 
of noise is known as measurement noise [ D| . 




Figure 1: Frequency spectrum of the exchange rate data y(t). For high 
frequencies the data are dominated by white noise. 

As will be seen in capter jO| , the existence of such an additive white noise 
dominating the small scales causes problems in our analysis, especially in the 
limit At — > which has to be performed according to equation (|7|). In order 
to avoid these problems, we applied a low-pass filter to the data. Each value 
y(t) was replaced by the weighted average of itself and its neighbouring data 
points, where the weighting function was chosen to be a gaussian centered 
at t with a width of 44sec. Figure |2] shows the result s(t) of this procedure 
in comparison to the original data y(t). 
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Figure 2: Samples of the time series of the original exchange rate data y(t) 
(cirles) and the smoothed signal s(t) (full line). 



According to equation fljjD , we can extract the noise n(t) by subtracting 
the smoothed signal from the original data y(t). If the assumptions leading 
to equation ([H|) are correct and the algorithm used to smooth the data is an 
appropriate one, the extracted n(t) should be white noise, i.e. it should be 
5-correlated with zero mean. 

The ratio of the mean value < n(t) > to the standard deviation a n of 
n(t) was found to be smaller than 8 • (cr n = 0.0003DM). This result 
justifies the assumption of a zero mean for n(t). The autocorrelation function 
Ru(t) of n(t) decreases from 1 to 0.08 within two seconds (see figure 
Compared to the temporal resolution of the data (the smallest time step 
between two subsequent data points is two seconds), R(r) indeed decreases 
rapidly. However, there appear to be small but significantly non-zero values 
for time delays r < 2mm, indicating that the approximation of the "real" 
signal s(t) by the smoothed signal is insufficient. We therefore restrict our 
analysis to time delays r larger than a certain elementary step T min which, 



from figure [3], we chose to be Amin. In units of r„ 
to be 5-correlated. 



R„ (t) can be considered 
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Figure 3: The autocorrelation function of the extracted additive noise n(t). 

Figure f| finally compares the probability density functions of price incre- 
ments x(t) calculated from the original data y(t) and the smoothed data s(t). 
It becomes evident that the influence of the noise n(t) is indeed restricted 
to small scales. For scales larger than r = 2r min the pdfs are practically 
identical. 

It may be worth to recall the result for fully developed turbulence, where 
a similar stochastic analysis was performed for velocity increments over space 
scales instead of price increments over time scales. In the case of turbulence, 
the Markovian properties were found to be valid only for (length-) scales 
and differences of scales larger than an elementary step size l mar [[UJ. But 
in contrast to the elementary step in turbulence which is a physical effect 
caused by the smoothing effects of viscous forces, the step r m j n used here for 
the analysis of financial data is due to the additional noise n(t) as discussed 
above. We want to mention that it is not unusual to find such a "cut off" 
at small scales for stochastic processes. For Brownian motion, the mean free 
path of molecules defines an elementary step size for the diffusion process in 
an analogous way. 

For the rest of the article, the price increments Xifa) are (unless it is 
mentioned explicitly) calculated from the smoothed data s(t). 
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Figure 4: Pdfs of price increments x(r) for r = 4mm, 8mm and 15mm (from 
bottom to top). The pdfs of x were calculated from the smoothed data (full 
symbols) as well as from the original time series y(t) (open symbols). Curves 
have been shifted in y-direction for clarity of presentation. Note that a ist not 
the respective standard deviation of each individual pdf, but the standard 
deviation of the whole data set y(t). 

3.2 Markov Properties 

Next, we give evidence that the Markovian property (§) holds. Strictly speak- 
ing, the relationships (Q) have to be verified for all positive values of N as 
well as for each set of scales T\, rjy, a task, which is evidently impossible. 
However, the data set considered here consisting of 10 6 samples allows to 
verify condition (^) for N — 3: 

p{x 1 ,Ti\x 2 ,T 2 ;x 3 ,T 3 ) = p(x 1 ,t 1 \x 2 ,t 2 ). (14) 

In figure ||, the contour plots of p(xi, T\\x 2 , r 2 ) and p(xi, r±\x 2 , r 2 ; x 3 = 0, r 3 ) 
have been superposed for the scales T\ = r min = 4min, t 2 = 2r min and 
r 3 = 3r min . The proximity of corresponding contour lines yields evidence 
for the validity of equation (|14D for the chosen set of scales. Additionally, 
two cuts through the conditional probability densities are provided for fixed 
values of x 2 . 
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Figure 5: (a): Contour plots of the conditional pdfs p(xi, Ti\x2, t 2 ) (dashed 
lines) and p(xx, t x \x 2 , t 2 ; x 3 = 0, r 3 ) (solid lines) for t x = r min , r 2 = 2r min and 

7~3 3T m j n . 

(b) and (c): Cuts through the conditional pdfs for x 2 = ±0.02cr. Open 
symbols: p(x\, Ti|x 2 , t 2 ; x 3 = 0, r 3 ), solid lines: Ti|x 2 , t 2 ). 

Figure |6] shows the same plots for a different set of scales: T\ = lh, 
T 2 = T i + T min and r 3 = t\ + 2r min . Again, we find a good agreement between 
corresponding contour lines. 
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Figure 6: (a): Contour plots of the conditional pdfs Ti|:r 2 , r 2 ) (dashed 
lines) and p(xx, Ti|x 2 , t 2 ; x 3 = 0,r 3 ) (solid lines) for T\ = lh, r 2 = T\ + r min 
and r 3 = ri + 2r min . 

(b) and (c): Cuts through the conditional pdfs for X2 = ±0.025a. Open 
symbols: p(x 1 ,t 1 \x 2 ,t 2 ; x 3 = 0,r 3 ), solid lines: p(xi, Tx\x 2 , r 2 ). 



Similar results were obtained for several other sets of scales chosen from 
the intervall r m j„ < T\ < 2h, with Ar = r 2 — ri = r 3 — r 2 ranging from r m j n 

to 3T m i n . 

Based on these result, we proceed with the well-founded assumption that 
the price increments of exchange rate data obey a Markov process for the 
range of scales under consideration, i.e. for scales and differences of scales 
larger than T min . 



3.3 Kramers-Moyal Coefficients 

According to equations (0) and @, the coefficients M k (x,r, At) can be 
calculated from the joint probability density functions. These joint pdfs 
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p(x, t — Ar ; x, r) are easily obtained from the data by counting the number 
N(x, x) of occurrences of the two increments x and x. 

According to equation (|7]), the limit Ar — > has to be performed next in 
order to obtain the Kramers- Moyal coefficients. Figure [j] shows the coefficient 
M 2 (x, r, Ar) for exemplarily chosen fixed values of x and r as a function of 
Ar. Throughout the intervall T min < Ar < 3r min , the dependence of M 2 on 
Ar turns out to be linear. For values of r < r min /2, the values deviate from 
that linear behaviour. It is interesting to note that the range of scales where 
those deviations begin is identical with the range of scales for which the 
autocorrelation function of the reconstructed additive noise n(t) has nonzero 
values (see figure |3|). The limit Ar — > is therefore performed by fitting a 
straight line to the M n in the intervall r m j„ < Ar < 2r min , thus avoiding 
problems with the values for Ar < r min /2 (see figure |7|). 

Figure [7] also displays the coefficient M 2 y \ which is obtained when instead 
of the smoothed signal s(t) the original data y{t) is used to calculate the 
coefficients M^. In the presence of additive white noise, the values of M 2 (y) 
diverge as Ar goes to zero. The limit Ar — > can not be performed in 
this case. Since the coefficients are nothing but conditional moments of 
the increments (see equation (|8|)), the reason for this behaviour can easily 
be understood by expressing the second moment of the increment of y(t) in 
terms of s(t) and n{t). Using equation (|T3|), we obtain: 

((y(t + r)-y{t)f) = ((s(t + r) + n(t + r) - s(t + r) - n(t + r) ) 2 ) 

= ((s(t + r)-s(t)) 2 ) + 2(n 2 ). (15) 

Due to the additive nature of the white noise n, the additional constant term 
2 ( n 2 ) arises which does not depend on the scale r. Similar terms also arise 
for the conditional moments Mfc when y(t) is used instead of s(t). Dividing 
the conditional moments of the increment by Ar according to equation (||) 
thus leads to the diverging behaviour of Mjf* in the limit Ar — » as shown 
in figure |7|. 

When the smoothed signal s(t) is used to calculate the coefficients M^, 
the linear dependence of M/,(x, r, Ar) on Ar is found to hold for k = 1,2, 
several scales r and all values of x. This enables us to calculate the coefficients 
Di(x, r) and D 2 (x, r) from the M\~ using the method of linear extrapolation 
described above. Figure |] shows the results for D\ and D 2 as a function of 
the price increment x at several scales r. Both coefficients exhibit simple 
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Figure 7: The coefficient M 2 (x, r, Ar) as a function of Ar for x = 0.05<r 
and r = 15mm (full circles), (a): logarithmic scale, (b): linear scale. For 
t > r mir) ,, M 2 is a linear function of Ar and can thus be extrapolated using 
a linear fit in the interval r m j„ < Ar < 2r m j n (full line). The open circles 
represent the results for M 2 for the same values of x and r when instead 
of the smoothed data s(t) the original data y(t) are used. The presence of 
additive white noise leads to diverging values for as Ar goes to zero. 



dependencies on the price increment. While Di is linear in x, D 2 can be 
approximated by a polynomial of degree two where the linear term is zero: 

Di(x,t) = -j(t)x, 

D 2 (x,r) = a(r) + (3(t)x 2 . (16) 
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Figure 8: The coefficients D\ (a) and D 2 (b) as functions of the price incre- 
ment x at scale r = 15mm (circles), r — 2h (squares) and r = 12h (crosses). 
The lines represent the results for D\ and D 2 obtained with the method 
described in section |3J^ (note the logarithmic presentation of D 2 ). 



Equation (jig) turns out to describe the dependencies of the coefficients 
Dk on x for all scales r up to two hours (see figure H). For larger scales, 
the statistics of the data are too poore to allow for a proper calculation of 
the coefficients M^, which results in considerable scatter of D\ and D 2 (as 
exemplarily shown for r = 12h in figure |8|). Indeed, from the analogous 
analysis of turbulent data fl6|| , we know that the conditional moments 



can be calculated properly only for scales r with r < 10 4 T, where T is 
the total length of the data set. With T being one year for the data set 
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considered here, we are thus restricted to scales smaller than one hour for a 
direct analysis of the Kramers- Moyal coefficients. 

Fitting the D k (x,r) by straight lines and parabolas, respecitvely, there- 
fore yields values for j(t), a(r) and (3{t) which scatter considerably and 
which do not exhibit a well defined functional dependence on the scale r. 
However, what can be concluded from the analysis of the coefficients is 
that the Kramers-Moyal coefficients are linear and quadratic functions of the 
price increment, respectively, with coefficients depending on the scale r. 

3.4 The Fourth Order Coefficient 

According to Pawula's theorem, it is of importance to estimate the fourth 
order coefficient and to decide whether it may be neglected. Figure || shows 
the coefficient M 4 (x, r, At) for fixed values of x and r as a function of At. 
Again, we find a linear dependence of the coefficient on At for At > T m i n . 
But whereas M 2 (Ar) increases as At goes to zero, M 4 decreases. Further- 
more, the linear extrapolation yields a value for which is small compared 
to the values of M 4 for T min < r < 2T min . This is a first hint that M 4 tends 
to zero in the limit At — > 0. 
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Figure 9: The coefficient M 4 (x,r, At) as a function of At for x = —0.05a 
and t = 15mm (circles). 

Performing the limit At — > for various values of x using the linear 
extrapolation confirms the above results. We obtain positive as well as neg- 



ative values for Z} 4 , as shown in figure |10] for the scale r = 15mm. Since, by 
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definition, D4 is a positive quantity, this means that the data set allows to 
take -D4 to be zero. 

2 10 6 1 , 1 1 1 

_ 110- 6 - 



E 




x / o 

Figure 10: The fourth coefficient D±(x, r) as a function of x for scale r = 
15mm (circles). The linear extrapolation yields positive as well as negative 
values for _D 4 . 

Similar results were obtained for several scales r < 2h. Based on these 
results, we will proceed with the assumption that the fourth order coefficient 
D4 vanishes, i.e. that Pawula's theorem applies and that the stochastic 
process is described by the Fokker-Planck equation. 

3.5 An Alternative Approach 

Based on the results of the preceeding section, we will proceed with the as- 
sumption that the statistics of the price increment is governed by the Fokker- 
Planck euqation with coefficients D\ and D 2 given by equation fll6|) and try 
to find estimates for the coefficients 7, a and (3. 

For the slope 7(1") of D\, such an estimate can be obtained from the 
conditional first order moment T x (r) = ( x(t) \xo(to) ) of the price increment. 
Multiplying the Fokker-Planck equation (|9]) with x, integrating with respect 
to x and finally using the result ([TBI) f° r D\, an equation for T x (r) can be 
derived. The result is: 



= -l{r)T\r). 



(17) 
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In deriving equation (|ITD , it was assumed that D 2 (x 1 r) p(x, t\xq, t ) goes to 
zero for \x\ — > oo. 

With the initial condition T l {r = r ) = (x(r = tq)\xo(tq)) = Xq, equation 
(|T7|) is solved by: 



T\r) 



Xq exp 



T(r') 



dr' 



;is) 



If the conditional moment T 1 follows the simple prediction (0), this can 
be taken as a strong hint for the validity of the assumptions underlying the 
derivation (validity of the Fokker-Planck equation and linear dependence of 
D\ on x). On the other hand, the simple dependence of T x (r) on 7 could 
also be used to measure the function j(t). 

Figure [11] shows the conditional expectation values T 1 /x as a function 
of t/t for scale r = 3h and several values of x . Indeed, as was expected 
from equation (p!8|), T 1 divided by xq turns out to be an universal function 
of t/tq. In particular, we do not find a systematic dependence on x$. 




t/t 



Figure 11: The conditional moment T 1 (divided by xq) as a function of the 
time scale r for r = 3h, xq = —0.05a (circles), xo = —0.025a (squares), xo = 

0. 05a (diamonds) and xq = 0.05a (triangles). The straight line represents a 
power law fit to the data yielding a scaling exponent of 0.95. 

Furthermore, it is found that T 1 is well described by a power law in r, 

1. e.: T 1 oc r c . It is easily seen that such a power law arises from equation 
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(18) for the case that j(t) is constant: 




T 1 (r) = xq exp \ / -j dr \ = Xq exp <j 7 In ( — 



(19) 

Hence, the conditional first order moment T 1 of the price increment x pro- 
vides a very convenient method to determine 7. Fitting the data presented 



in figure [LI], we obtain a value of 7 = 0.95 for the scaling exponent. Summa- 



rizing the results obtained by varying the values of xq and r finally yields: 

7 = 0.93 ±0.02 . (20) 

It is interesting to note that the value of 7 we obtain is close to one. However, 
the statistics of the data set we used is rather poor and further investiga- 
tions using larger data sets will be necessary to decide whether 7 deviates 
significantly from one. 

Estimates for the coefficients a(r) and (3{t) of D2 can be obtained using 
the equations for the n-th order moments < x{r) n > of the price increment. 
With the result ( filf ) for D\ and D 2 , equation flT2]) yields for n > 2: 



-rj-(x n (r)) = - ni (x n (T)} + n(n-l)a(T) (x n - 2 (T)) 



dr 



which can be rewritten as: 



+ n(n-l)/3(r) (x n (r)) , (21) 



T 



9 



n(xHr)) y-^-mr)-{n-l)a(r)^^. (22) 

As 7 is known from the analysis of the conditional first order moment, two of 
the equations ( p2"D are sufficient to calculate the two unknown functions u(t) 
and /3(t) from the moments (x n (r)). We use the moments of order two and 
four, since even moments of low order can be determined best from empirical 
data. Due to the rather limited number of samples it is even uncertain 
whether the moment of order six can be determined properly (see [IE] for a 
more detailed discussion of this topic). 
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Figure (|T^) shows the moments (x n (r)) for order two and four as functions 
of the time scale r. Throughout the range r m j n < r < 12h, the moments 
can be described by power laws in r: 

(x n {r)) oc r c ". (23) 

Using this result, the left hand side of equation (p2|) can easily be shown 
to be equal to — . Fitting the empirically determined moments with power 
laws, we find the following values for the scaling exponents ( n : 

( 2 = 0.94 ± 0.04 

C 4 = 1.72 ± 0.12. (24) 

The estimates for the errors were obtained by varying the range of scales used 
for the fit. Note that we found a value of (2 which is slightly smaller than 
one, the value which is usually assumed PUL M. Due to the poor statistics of 
the data set, however, this result is by no means significant; the second order 
moment can be fitted with a scaling exponent of (2 — 1 with almost the same 
accuracy. Yet, the values given in (|24]) are the best fitting coefficients for the 
data set under consideration and will therefore be used here. 

10- 1 
10- 3 

A 
c 

2 10 4 

X 

V 10- 

10 s 

10- 7 

t [sec] 

Figure 12: The moments (x n (r)) of the price increment x(r) as functions of 
the time scale r for order n = 2 (top) and n = 4 (bottom). Throughout the 
range r min < t < 12h the moments are in good approximation described 
by power laws (straight lines). 
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Using the results for the second and fourth order moments, equations (|22|) 
can easily be solved and yield estimates for the unkown functions u(t) and 
/3(t). We find that a is a linear function of the time scale r with a slope «o 
of 0.016d~ l , while f3 is approximately constant in r and has a value of 0.11. 
Let us summarize our results for Dx and D 2 : 

Di(x,t) = -71 

7 = 0.93 ± 0.02 

D 2 (x,r) = a r + (3x 2 

a = 0.016CT 1 ± 0.002CT 1 

P{t) « const = 0.11 ± 0.02. (25) 



3.6 Consistency Checks 

The results and assumptions of the preceeding sections lead to Fokker- Planck 
equations for the pdf p(x,r) and the conditional pdf p(x, t\xq, to), respec- 
tively (equations ( |T0D and (§)). The coefficients Di and D 2 which completely 
determine these equations were estimated from the data. Thus we can claim 
to have found the complete description of the stochastic process for the data. 
In order to test this result, we compare the (numerical) solution of the Fokker- 
Planck equation with the distributions obtained directly from the data. The 
algorithm used for the numerical iteration is based on the approximative so- 
lution of the Fokker-Planck equation for small steps At. According to [14 . 
the conditional pdf p{x\, t — At\xq, t ) is, for small At and arbitrary Di and 
D 2 , a Gaussian distribution in x\ with mean value xq — D\At and standard 
deviation \J2D 2 At. In order to obtain the conditional densities for larger 
steps, we use the Chapman-Kolmogorov equation 

+00 

p(x 2 ,t -2At\x ,t ) = J p(x 2 ,T -2AT\x 1 ,T -AT)p(x 1 ,T -AT\x ,T )dx 1 . 

— CO 

(26) 

( |26| ) is a direct consequence of the Markov condition (f|). Iterating this proce- 
dure, we finally obtain the conditional pdf p(x, r — uAt\xq, t ). Multiplying 
with p(xq, t ) and integrating with respect to xo yields the pdf p(x, t — uAt). 

Figure [13], which compares the solutions of the Fokker-Planck equation for 
the pdf p(x, t) with the empirically estimated pdfs, proves that the Fokker- 
Planck equation accurately describes the evolution of p(x, r) in r over the 
range 4min < r < 12h. 
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Figure 13: Comparison of the numerical solutions of the Fokker-Planck 
equation (solid lines) for the pdfs p(x, r) with the pdfs obtained directly 
from the data (open symbols). The scales r are (from top to bottom): 
r = 12h, 4h, lh, 15mm and 4mm. The pdf at the largest scale r = 12h 
was parametrized (dashed line) and used as inital condition for the iteration 
of the Fokker-Planck equation. Curves are shifted in vertical direction for 
clarity of presentation. 

As mentioned above, the Fokker-Planck equation also governs the condi- 
tional pdf p(x,t\xo,t ). As a further test of our results, we calculated the 
solutions of the Fokker-Planck equation (|9|). Figure [14] shows the result for 
ro = lh and r = 0.5/i, again in comparison with empirical data. Taking into 
account the various uncertainties and assumptions in the determination of 
the coefficients D\ and D2, the agreement between the solution of @ and 
the data is remarkably good. 
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Figure 14: Comparison of the numerical solution of the Fokker-Planck equa- 
tion for the conditional pdf p(x, t\xq, tq) with the emprical data. 

(a) : Contour plots of p(x, t\xq, tq) for ro = lh and r = 0.5h. Dashed lines: 
numerical solution of eq. (|]), solid lines: empricial data. 

(b) and (c): Cuts through p(x, r\x , r ) for x = — 0.062cx and x = +0.062cr, 
respectively. Open symbols: empirical data, solid lines: numerical solution 
of the Fokker-Planck equation. 



4 Summary and Conclusions 

The purpose of the present article is to show how the mathematical frame- 
work of Markov processes can be applied to the analysis of emprical high 
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frequency exchange rate data. Within the limitations due to the finite num- 
ber of samples, we were able to verify the Markovian properties of the price 
increment for time delays r larger than r m j n = 4min. 

Furthermore, we obtained estimates for the Kramers-Moyal coefficients 
-Dfc(x,r) and found hints that the fourth order coefficient D 4 is zero. The 
comparison of the solutions of the resulting Fokker- Planck equations with 
the empricial distributions strongly supports our results for D\ and D 2 . 

It is worth to note that the results given in the present article differ from 
the results of our previous analysis presented in JL2]|. While in (12| D\ and 



D 2 were also found to be linear and quadratic functions of x, respectively, 
the values given for the coefficients 7, a and ft were different from the ones 
given here in equation 



The analysis presented in is incomplete, because it did not take into 
account the presence of the additive white noise (see section |3.1| ) and mainly 
focused on the functional dependence of the on x. The numerical values 
of the coefficients 7, a and ft were chosen such that the pdfs of the price 
increment x were fitted best by the resulting Fokker-Planck equation. How- 



ever, it can be shown that the coefficients given in [12 1 fail to describe the 
conditional pdf p(x, t\xq, t ), i.e. cannot be correct. 

This is a quite important result: There are several choices for the coef- 
ficients which correctly describe the one-point distribution p(x, r), but only 
one choice which also describes the joint statistical properties of several in- 
crements on different scales. This, in turn, means that a complete statistical 
characterization of the price increments (or returns) necessarily has to in- 
clude the multi-scale statistics and must not be restricted to the analysis of 
increments (or returns) on only one scale. 

The results of our analysis presented in this paper also eliminate the 
discrepancy between our approach as presented in [T2] and the often reported 



power-law behaviour of the distributions of returns. As previously shown in 
I?!, the solutions of the Fokker-Planck equation asymptotically show scaling 



behaviour with a scaling exponent /i which can be calculated from the linear 
and quadratic coefficient of Di and D 2 , respectively: 

" = W (27) 



While the values reported in |L2] lead to a value of about 12 for the scaling 
exponent /1, in disagreement with the usually observed value between 3 and 
5, our new values (|25D yield /i = 4.2 ± 0.8. 
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Another interesting result of our analysis concerns the symmetry of the 
stochastic process underlying the evolution of the price increment x in the 
scale variable r. According to (|25|), the drift coefficient D\ is completely 
antisymmetric in x (Di(—x) = —Di(x)), while the diffusion coefficient D 2 is 
symmetric in x. The resulting Fokker-Planck equation fllQ|) for the pdf p(x, r) 
is symmetric in x. This means that if the large scale pdf, which serves as the 
initial condition for the partial differential equation (|H]), was symmetric, it 
would remain symmetric for all scales r. In other words: The asymmetries 
of the distributions of price increments on small scales are solely due to the 
asymmetry of the increment at the largest scale, i.e. they are a consequence 
of some market event on long terms. The stochastic process itself, at least 
on scales smaller than a day, is symmetric in x. This result might provide a 
possibility to distinct between trends (as for example caused by the long term 
development of a country's economy) and the short term fluctuations caused 
by the market itself i.e. by the interactions of several agents trading an asset. 
This fits well with the results of commonly discussed market models, which 
are quite succesfull in qualitatively explaining most of the features of real 
markets but which do not naturally exhibit asymmetries between positive 
and negative returns, see for example and JHJ . Our result indicates that 
there is no need to incorporate asymmetries in such models. However, so 
far we can only present results for a single data set over an one-year period. 
Further studies on a larger variety of assets and time scales will be necessary 
to check the significance of this particular result. 

Finally, we would like to comment on the proposed analogy between fi- 
nancial markets and turbulent hydrodynamic flows, cf. [p0| . The complex 
statistical behaviour of velocity increments on a certain length scale in turbu- 
lent flows is assumed to be due to a cascading process. Energy, which is fed 
into the system on large scales, is continuously transported towards smaller 
scales due to the inherent instability of vortices of a given scale towards per- 
turbations on smaller scales. Finally, the energy is dissipated at the smallest 
scale. 

A similar mechanism has been proposed for financial markets, where the 
energy cascade was replaced by a flow of information. Initially, the assump- 
tion of a cascading process in financial markets was based on similarities in 
the empirical description of the pdfs of price and velocity increments, respec- 
tively [pOR . The analogy between turbulence and finance has inspired many 



further studies (see, for example, [^1|], |fffi| , [p3| ) ; but has also been criticized 
for being too superficial M. 
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However, the results of the Markov analysis of financial and turbulent 



data (see [16] for details on the Markov analysis of turbulence) allow for 
a more detailed comparison of the stochastic processes underlying the two 
phenomena with the predictions of common cascade models. 

The most prominent model for turbulence is Kolmogorov's theory dating 
back to 1962 (furtheron refered to as K62). The prediction of the K62 model 
for the coefficients Di and D 2 can easily be shown to be |24| : 

D 1 = -jx, 

D 2 = (3x 2 , (28) 



where 7 and (5 both are constant (7 ~ 1/3 for turbulence). Castaing 



proposed the model of a multiplicative cascade where the increment on a 
given scale is linked to an increment at a larger scale via a random multiplier. 



This model was shown to imply Markovian properties |26| and to predict 
coefficients D\ and D 2 similar to those predicted by K62, except that it 
allows for a scale dependence of 7 and (3. 

Recently, Dubrulle []27| proposed an extension of Castaing's multiplicative 



cascade model which incorporates additive noise. Comparing the results (p5|) 
of our analysis with the predictions of the various models, we can state that 
the statistics of the price increment is well described by the model of a 
multiplicative cascade with additive noise (or, in the notion of Jl7|], by the 
mechanism of "multiplicative noise with reinjection"). 

It is interesting to note the result of the Markov analysis for fully de- 
veloped turbulent flows JIB]] . It turns out that measured turbulent data are 
not completely described by any of the above mentioned models. While D\ 
turns out to be linear, D 2 exhibits are more complicated dependence on the 
increment than predicted by the multiplicative cascade models. In partic- 
ular, we find a linear term in D 2 which breaks the symmetry and which is 
by no means negligable, see IIBJ. In some sense, one might therefore say 



that financial data are much closer to ideal " Kolmogorov-like" turbulence 
than turbulent data. Besides this analogy we found a clear differene in the 
mechanisms leading to the finite step-size of the cascade processes. In the 
case of turbulence a smoothing out due to viscosty was found, while financial 
markets seem to be dominated by some uncoorelated additional noise on the 
smallest scales. 

Summarizing, it is the concept of a cascade in time hierarchy that allowed 
us to derive the results of the present paper, which in turn quantitatively 
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supports the initial concept of an analogy between turbulence and financial 
data. Furthermore, we have shown that the smooth evolution of the pdfs 
down along the cascade towards smaller time delays is caused by a Markov 
process with multiplicative noise. 

Metaphorically, our characterization of volatility of currency markets as 
random and uneven flow of information to the markets is similar to |7| which 
describes stochastic process of the stochastic volatility model as random and 
uneven flow of news to the financial markets. 

Acknowledgements 

We gratefully acknowledge useful discussions with A. Soofi, D. Sornette, Zhi- 
Feng Huang, St. Luck and M. Siefert. The FX data set has been provided 
by Olsen & Associates (Zurich). 

References 

[1] L. Bachelier, Theorie de la speculation, Ph.D. Thesis, Annales Scien- 
tifiques de l'Ecole Normal Superieure III- 17, 1900. 

[2] Engel, R. F. and M. Rothschild, (eds.). (1992). Journal of Econometrics, 
A supplement to the Journal of Econometrics, Vol. 52. 

[3] Bollerslev, T., 1986. Generalized autoregressive conditional heteroskedas- 
ticity. Journal of Econometrics. 31, 307-327. 

[4] Bollerslev, T., Chou, R. Y., Kroner, K.F., 1992. ARCH modeling in fi- 
nance: a review of the theory and empirical evidence. Journal of Econom- 
terics, 52, 5-59. 

[5] Taylor, S. J., 1994. Modelling stochastic volatility: a review and compar- 
ative study. Mathematical Finance, 4, 183-204. 

[6] Ghysels, E., Harvey, A., and Renault, E. (1996). Stochastic Volatility. 
In: Handbook of Statistics (vol. 14), eds. G. S. Maddala and C.R. Rao, 
Amsterdam: North-Holland, pp. 119-191. 

[7] Gallant, A. R., Hsieh, D., Tauchen, G., Estimation of stochastic volatility 
models with diagnostics. Journal of Econometrics, 81, 159-192. 



Renner, Peinke, Friedrich 



28 



[8] R. Mantegna & H.E. Stanley 1996. Turbulence and financial markets. 
Nature 383, 767. 

[9] R. Mantegna & H.E. Stanley 1999. An introduction to econophysics. 
Cambridge University Press, Cambridge. 

[10] H.E Stanley et al., Quantifying fluctuations in economic systems by 
adapting methods of statistical physics, Physica A 287 (2000) 339 - 361. 

[11] Voit, J., The Statistical Mechanics of Financial Markets. Springer- 
Verlag, Heidelberg (2001). 

[12] Friedrich, R., Peinke, J. and Renner, Ch. (2000). How to quantify de- 
terministic and random influences on the statistics of the foreign exchange 
market. Phys. Rev. Lett, 84, 5224-5227. 

[13] Heald, J.P.M. and Stark, J. (2000) Estimation of Noise Levels for Models 
of Chaotic Dynamical Systems. Phys. Rev. Lett., 84, 2366-2369 

[14] Risken, H., The Fokker-Planck equation, ( Springer- Ver lag Berlin, 1984); 
Van Kampen,N. C, Stochastic processes in physics and chemistry (North 
Holland, Amsterdam, 1981); Gardiner, C.W. Handbook of Stochastic 
Methods, ( Springer- Verlag Berlin, 1983) 

[15] Hanggi, P. & Thomas, H. 1982 Stochastic processes: Time evolu- 
tion, symmetries and linear response Physics Reports 88, 207. 

[16] Rener, Ch., Peinke, J. and Friedrich, R. Experimental indications for 
Markov properties of small scale turbulence. JFM 433 (2001), 383-409. 

[17] D. Sornette, Fokker-Planck equation of distributions of financial returns 
and power laws. Physica A 290 (2001) 211-217. 

[18] R. Cont, J. P. Bouchaud, Herd behaviour and aggregate fluctuations in 
finacial markets. |cond-mat /9712318 . 

[19] D. Stauffer, D. Sornette, Self-organized percolation model for stock mar- 
ket fluctuations, Physica A 271 (1999) 496-506. 



[20] Ghashghaie, S., Breymann, W., Peinke, J., Talkner, P., Dodge, Y., 1996. 
Turbulent cascade in foreign exchange markets. Nature 381, 767-770. 



Renner, Peinke, Friedrich 



29 



[21] A. Arneodo, J.-F. Muzy & D. Sornette, "Direct" causal cascade in the 
stock market. Eur. Phys. J. B 2 (1998), 277-282. 

[22] F. Schmitt, D. Schertzer & S. Lovejoy, Multifractal analysis of foreign 
exchange data. Appl. Stoch. Model D. A. 15 (1), 1999. 

[23] W. Breymann, S. Gashghaie & P. Talkner, A stochastic cascade model 
for FX dynamics. International journal of theoretical and applied finance 
(IJTAf) 3 (3), 357-360, 2000. 

[24] Friedrich, R., Peinke, J., Physica D 102, 147 (1997); Friedrich, R., 
Peinke, J., Phys. Rev. Lett. 78, 863 (1997) 

[25] B. Castaing, Y. Gagne & E.J. Hopfinger, Physica D 46, 177-200 (1990); 
B. Castaing, J. Phys. II France 6 105 (1996). 

[26] P.-O. Amblard & J.-M. Brossier, On the cascade in fully developed tur- 
bulence. The propagator approach versus the Markovian description. Eur. 
Phys. J. B 12 (1999), 579-582. 

[27] B. Dubrulle, Affine turbulence. Eur. Phys. J. B 13, 1 (2000). 



